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ON A NON-HOMOGENEOUS AND NON-LINEAR HEAT 

EQUATION. 

LUCA BISCONTI AND MATTEO FRANCA 


Abstract. We consider the Cauchy-problem for a parabolic equation of the 
following type: 

^ = Au + f{u, |a:|), 

where / = f(u,\x\) is supercritical. We supply this equation by the initial 
condition u(x, 0) = (p, and we allow (p to be either bounded or unbounded in the 
origin but smaller than stationary singular solutions. We discuss local existence 
and long time behaviour for the solutions u(t,x](p) for a wide class of non- 
homogeneous non-linearities /. We show that in the supercritical case, Ground 
States with slow decay lie on the threshold between blowing up initial data 
and the basin of attraction of the null solution. Our results extend previous 
ones allowing Matukuma-type potential and more generic dependence on u. 

Then, we further explore such a threshold in the subcritical case too. We 
find two families of initial data C(^) 'ip{x) which are respectively above and 
below the threshold, and have arbitrarily small distance in L°° norm, whose 
existence is new even for f{u,r) = . Quite surprisingly both C{x) and 

'ip{x) have fast decay (i.e. while the expected critical asymptotic 

behavior is slow decay (i.e. 


Key Words: Cauchy-problem, semi-linear heat equation, singular solutions, 
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1. Introduction 

The purpose of this paper is to study the asymptotic behavior of positive solu¬ 
tions of the following Cauchy problem 

(1-1) ^ +f{u,\x\), 

(1.2) u{x,0) = 

where x G R", n > 2, and / = f{u, |x|) is a potential which is null for u = 0. 

In the last 20 years this problem has raised a great interest, starting from the 
model cases f{u, |a;|) = and f(u, lx]) = \x\^u'^~^. Due to symmetry reasons, 
from now on we use notations which are standard for the stationary problem, so 
we refer to /(u, |a;|) = as the model case. In so doing it will be more clear 
the relationship between the critical values for 113 appearing below, and their 
meaning in other contexts of functional analysis. 

We assume that / is supercritical with respect to the Serrin critical exponent, 
i.e. 2* := , and for some specific results we require / to be supercritical also 

to the Sobolev critical exponent, i.e. 2* := The exponents 2* and 2* are 

related to the continuity of the trace operator in L'^ and to the possibility to embed 
in , respectively. 

Here, we want to analyze the structure of the border of the basin of attraction 
to the null solutions, and the set of initial data </) of solutions of (inD-diU) which 
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blow up in finite time. Our main aim is to extend the discussion to a wide class of 
potentials: For the remainder of the paper we will always assume the following 
FO: The function /(u, r) is locally Lipschitz in u and r for any u > 0 and 
r > 0. Moreover /(0,r) = 0, f{u,r) > 0 and f{u,r) is increasing in u, 
for any u > 0 and any r > 0, and there is a constant C{u) > 0 such that 
f{u,r)r‘^ < C{u) for 0 < r < 1. 

Further hypotheses on / will be given in the sequel (see conditions GO, Git, and 
Gs in Section (2). Possible examples are the following 

(1.3) /(m, |a:|) =/ci(|x|)|u|«i"\ 

(1.4) /(m, |a:|) = ki{\x\)\u\'^^~^ + fc 2 (|a;|)|u|® 

(1.5) f{u, |a:|) = fci(|a;|)min{u'Ji“\u'’=“^}, 

where qi < q 2 and ki = ki{\x\), i = 1,2, are supposed non-negative and Lipschitz 
continuous, and such that 

(1.6) ki{r) ~ AiT^' as r —0, ki{r) ~ as r ^ -foo, 

where Ai, Bi > 0, -^i > 0, > 0, > 2 and 5i, rji > {n - 2)q - {n- 1), for 

i = 1,2 (so for 6i = rji = 0 we require qi > 2,). More precise requirements on ki, 
i = 1,2, will be provided later on according to the ones on /. 

Due to the nature of the considered potentials, in general we cannot expect the 
solutions of to be differentiable, or even continuous, everywhere. In fact, 

we deal also with solutions that may be not defined at a; = 0 since they become 
unbounded. 

In Section [3] we prove the existence of a proper class of weak solutions to the 
considered problem (see Lemma 13.71 and Theorem 13.81 below 1 and we actually show 
their improved properties. We consider the classes of Gs-mild and Gg-mild solu¬ 
tions to (ini)-(irii) ( see the definitions 1321 and 1321 below, see also m) proving 
local and global existence as well as uniqueness. 

Let u{x, t; (j>) be the solution of (ll.ll) - (ll.2l) . The analysis of the long time behavior 
of u{x, t] (j)) is strongly based on the separation properties of the stationary solutions 
of (II.m . i.e. functions u{x) solving 

(1.7) Au + f{u,\x\) = 0, 

and in particular on the properties of radial solutions. Notice that if u{x) is a radial 
solutions of ini), setting U{r) = u{x), for r = |a;|, then U = U{r) solves 

77 — 1 

(1.8) [/" +- U' + f{U,r) = 0, 

r 

where “' ” denotes the derivative with respect to r. 

In the whole paper we use the following notation: U (r) is regular if [/(O) = a > 0, 
so we set t/(r) = U{r, a), and we say that t/(r) has a non-removable singularity (or 
shortly that it is singular) if limr_>.oC/(r) = -too. Similarly, we say that a positive 
solution V{r) of p.8ll has fast decay (f.d.) if hmf._>.+ooI^(»")?’"'~^ = /3 > 0 and we set 
V{r) = V{r,j3), and that V{r) has slow decay (s.d.) if hmj._>+ooF(r)r"'“^ = -|-oo. 

Further, U{r) is a ground state (G.S.) if it is a regular solution of (11.81) which 
is positive for any r > 0. Instead, we say that U{r) is a singular ground state 
(S.G.S.) if it is a singular solution of (|1.8D which is positive for any r > 0. The 
asymptotic behavior of singular and slow decay solution is well understood and will 
be discussed in more details in Section [2] 

Roughly speaking, the w-limit set of dm is (usually) made up by the union of 
solutions of Equation (II. 8L see e.g. [laiiniiii], and these solutions are one of the 
ingredient to construct sub and super-solutions to (HD , see e.g. 
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We briefly review some known results concerning ini) and We need 

to introduce some additional parameters which play a key role in what follows: 
Recalling that 2* = 2^^ 


and that 2* = we have 


(1.9) 


Pf < 2* < 2* < cr* where Pp ■= 2 


n + 1 


a 


* 


(n—2)^—4n+8\/n—1 
(n-2)(n-10) 
+00 


if n > 10 
if n < 10 


The parameters Pp < 2* < 2* are critical exponents for this problem and their role 
will be specified few lines below. Here, Pp — 1 is the so called Fujita exponent. 

We assume first that / is of type la, with fci(|a;|) = ko{r), r = |a;|, and 
ko{r) := Kor^°, where Kq is a positive constant and > —2. Let us also introduce 
the followings 


( 1 . 10 ) 


lo 2 


g + ^0 

2 + Sq 


and 



2 Sq 
g-2’ 


and mo = mo(Zo)- 

In this case, whenever Iq > 2*, we have at least a S.G.S. with s.d. (j)s{x) = 
where Pi > 0 is a computable constant, which is unique for Iq ^ 2*. 
Also note that if Jq = 0 then Iq = q and ttiq = 2j{q— 2). Moreover, all the regular 
solutions of (II. 8|) have a non-degenerate zero for Iq £ (2, 2*), they are G.S. with f.d. 
for /o = 2*, and they are G.S. with s.d. for Ig > 2* (see, e.g. [H]). 

Again, if 2* < Ig < a*, then all the regular solutions cross each other, while if 
Ig > a* and a 2 > ai, then U{r,a 2 ) > U(r,ai) for any r > 0, see [57]. In fact, 
when the structure of positive solution of m changes, the asymptotic behavior 
of solutions of (ll.ll) - (ll.2l) changes too. 

Let us recall that all the solutions of (fTH-lfTTl) blow up in finite time if Ig < Pp, 
so the null solution is unstable in any reasonable sense (see [B IlSj). If Ig > Pp 
the null solution is stable with the suitable weighted L°°-norm, but still “large” 
solutions blow up in finite time. 

There are several papers devoted to explore the threshold between the basin of 
attraction of the null solution and the set of initial data which blow up in finite 
time (see, e.g. |571|T5J|T3] andalso [1911501 [5T| 1. It seems that radial G.S. of dni) 
play a key role in defining such a border. In particular Gui et al. in m (see also 
[27]) proved the following: 


1.1. Theorem. [T5] Assume f{u,r) = Kgr^°u‘‘~^ and 2* < Ig < a*. Then 

(1) : If 4>{x) ^ U(\x\,a) for some a > 0, then ||M(a:, t; ^(')||oo —>-0 as t ^ -l-oo. 

(2) : If 4>{x) ^ U{\x\,a) for some a > 0, then u(x,t;(j)) must blow up in finite 

time, i.e. there is T^j, £ (0, oo) sueh that limt-^.j',^ ||'w(x, i;</')l|oo = -1-00. 


This result was extended in [I] to potentials / of the form where 
(1.11) ki{\x\) = kg{r) := r^°K{r) 

with K{r) varying monotonically between two positive constants. Then, in [28] it 
was extended to / of the form (11.41) where ki{r) = r^'ki and fci > 0 is a constant. 
An interesting related topic is the rate of decay of fading solutions and of blow up 
(see, e.g. [571I5]). 

It is worth mentioning that in Ezi [ig Wang et al. proved that for f{u) = 
the situation is very different if g > cr*, i.e. G.S. are stable and weakly 
asymptotically stable with the suitable weighted L'^-norm. This result still holds 
true also for / as in (O if kg{r) = r^K{r), when K is decreasing, uniformly 
positive and bounded, and lg> a*, where Ig is defined in (11.101) . The same result 
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holds for (HI) when fci = 1, for i = 1, 2 and q 2 > Qi > cr* (see [IH])- The extension 
of this stability results to the potentials considered in this paper will be object of 
future investigations. 

A first contribution of the present paper is the extension of Theorem 11.11 to a 
number of non-linearities including da and a- In fact, we propose a unifying 
approach that allows us to consider a wider class of non-linearities including, e.g. 
a among others. 

As we have already seen, the sub- and supercriticality of (11.11) in the non- 
homogenous case, e.g. if the potential / is as in a, depends on the interplay 
between the exponent q = qi and the asymptotic behaviour oi k = ki. The same 
happens for the asymptotic behavior of positive solutions to (11.81) . Therefore we 
define the following parameter, useful to combine the two effects: 

(1.12) I = l{q^5)-=2^^^ and m{l) := 

If f{u,r) = as in the Wang case, we have a subcritical behavior for I < 2* 

and supercritical behavior for I > 2*] the same happens with the other critical 
parameters defined in (11.91) . 

We stress that singular and slow decay solutions U (r) of Equation (11.81) behave 
as ~ as r —>■ 0 and as r —)► -|-oo respectively (where Pi is a computable 

constant). 

Using different values of I we can allow two different behaviors for singular and 
slow decay solutions, namely: Denote by lu and m(lu) the parameters ruling the 
asymptotic behavior of singular solutions ?7(r), i.e. U{r) ^ g^g ^ q- 

similarly, Ig and m{ls) are the parameters ruling the asymptotic behavior of slow 
decay solutions V{r), i.e. V{r) ~ as r —>• -|-oo. 

This will allow us also to consider Matukuma potentials (see below, and see also 
[25]1: Thus, e.g. in the case dESD with fci = fc as in (HH), we have 

{ / — ( 1 = 0 9+’? 

milu) = ^ [ m{lg) = ^ 

Analogously, in the cases (11.41) and (11.51) we have, respectively, that 

/ as in (fUdl) . then l„ = max{2^^ | i = 1,2}, 4 = min {2^^ | i = 1,2} 
^ ^ / as in (H^, then 2^^, 4=2^^ 

and according to (I1.12|) we also obtain m{lu) = and to(U) = 

Let us state the following sub and super-criticality conditions related to fc^r), 
1 = 1,2, that replace the fact that, for fco(r) = r^°K{r) in (11.111) . K{r) is monotone: 

H+ Jo > 0 for any r > 0 and any i, strictly for 

some i and r > 0. 

■ /o^< 0 for any r > 0 and any i, strictly for 
some i and r > 0. 

We emphasize that if / is either of type (11.81) . (11.41) or (11.51) if iT-|- holds then 
regular solutions of (11.81) are crossing while, if H— is verified, then they are G.S 
with slow decay (see 0). 

Now we can state the following: 

1.2. Proposition. Assume that / is either of the form mw and 

satisfies H—. Further, assume > 2*, and 2* < Ig < a*. Then all the regular 
solutions U{r,a) of il.8\} are G.S. with s.d., and there is at least a S.G.S. with 
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slow decay U{r,oo). Moreover if 0 < ai < a 2 < oo for any M > 0 there is R = 
R{a 2 ,ai) > M such that U{R,a 2 ) = U{R,ai) and -^U{R,a 2 ) — -^U^R, ai) < 0. 

This result is a direct consequence of Proposition 12.121 and Remark 12.101 below. 
In fact, the intersection property of G.S. is a secondary contribution of this paper. 

In this setting we can extend Theorem 11.11 as follows 

1.3. Theorem. Assume that the assumptions of Provosition are verified, then 
the same conclusion as in Theorem \l.l\ still holds true. 

The above result is obtained as a corollary of Theorem 14.11 below, which is 
somewhat more general. 

We highlight the fact that when / is of type Theorem 11.31 generalizes the 
result of [I] to the case where fci(r) is not monotone decreasing and may even be 
increasing in some cases. E.g., let / be of type (IE3D with k{r) = fci(r) = 1 + r“; 
assume q > 2* and a > 2*(g — 2) — 2, so that from (11.121) we have lu = q and 
Is = 2{q + a)/{2 + a) > 2*, then Theorem 11.31 applies directly to this situation. 

Notice that Theorem 11.31 requires a weaker condition on /„ than on Ig. Hence, 
Theorem 11.31 applies also to the case (11.31) even for q> a*, with the condition that 
a € ( ,^_^iq cr*), ” 2^(9 2*)], while from [57] and n we know that in this case, 

if k(r) is a constant or a decreasing function varying between two positive values, 
G.S. are stable, so we are in the opposite situation. 

Also, we emphasize that Theorem 11.31 extends [TJ Theorem 1] also to Matukuma 
type potential (see, e.g. ESI for more details), which are a model in astrophysics, 
i.e. to / of the form (|1.3I) where q £ [2*,cr*) and k{r) = 1/(1 + r“), where a G 
{0,2-a*{q-2)). 

When / is of type (11.41) we extend the result in [55] to the case where ki{r), 
1 = 1,2 are r-dependent functions, and we can deal with a generic family of non- 
linearities including dEll). 

Let us go back again to the case of f{u,r) = f{u) = The singular solu¬ 

tion 4>s{x) := Pi\x \~'^/seems to play a key role in determining the threshold 
between solutions converging to zero and solutions blowing up in finite time. 

In [13 Ni shows that if 2* < g < cr* and (j){x) < 4>s{x), then u{x,t;cj)) converges 
to the null solution as t —>■ -boo. Let Ai denote the first eigenvalue of the Laplace 
operator in the ball of radius r = 1; if liminf| 2 ,|_^_|_oo cj){x)\x\~'^^^‘>~‘^^ > 
then u{x, t; (f) blows up in finite time. 

Wang in m shows that if g > cr* and liminf| 2 ;|_> +oo(/)(x)|x| > Pi then 

u{x,t;(f>) blows up in finite time. Note that this result is optimal since, for g > 
cr*, there are uncountably many G.S. with s.d. asymptotic to Pias 
|a;| —> -boo. On the other hand, in EZl Theorem 0.2, point (ii)], Wang proved the 
following. 

1.4. Theorem. Consider f{u,r) = f{u) = u'^~^, where 2* < g < a*; then for any 

P > 0 there is a radial decreasing upper solution of (13) Xgix) such that X/ 3 ( 0 ) = /3 
and —>■ Pi ;= [m(g)(n — 2 — TO(g))]^/('^“^^ as |a;| —>■ -boo. Moreover 

limt_,.+oo \\u{x,t-,xp){^ + |a;r)||oo = 0 for any 0 <n < m{q). 

In fact, the result is proved for a slightly more general potential f{u, r) = r^u'^~^, 
with a proper S. These results seems to indicate (or more in general the 

decay rate of slow decay solutions of (LTD) as the optimal decay rate for having 
solutions which are continuable for any t £ M., see the introduction of [13] for a 
detailed discussion on such a topic. 

From now till the end of this section we consider / as follows: 

(i): / as in m and ki satisfies m with lu,ls > 2*; 
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(ii): / as in dEl and ki, i = 1,2, satisfy (fra with lu,ls > 2^ ; 

(iii): / as in (11.51) and k satisfies (HID with lu,ls > 2*; 

As a a consequence of Theorem 14.21 we generalize this result to present setting, i.e. 

1.5. Theorem. Assume f either of the form (i), (ii), or (iii). Assume either H-\- 

with lu,ls G (2*, 2*] or H— with lu ^ , Is & [2*,cr*). Then we have the same 

conclusion as in Theorem \1.4\ but m[q) is replaced by m{ls) and Pi is replaced by 
the computable constant (e.g. P^°° = [m{ls){n — 2 — iff is of 

type (i);. 

We emphasize that, as far as we are aware, this result is new anytime we consider 
f{u,r) as in (|1.3D but fc(r) ^ r^, and for (11.41) even for ki = k 2 = 1. 

The main contribution of this paper is the following result (consequence of the 
slightly more general Theorem l4.dll which goes in the opposite direction with respect 
to Theorem 1 1. 41 ('and hence to Theorem ll.511 . and shows that the situation is really 
delicate. 

1.6. Theorem. Assume f either of the form (i), (ii), or (iii). Further assume that 

either lu,ls O'f'e. in (2*, 2*], and H-\- holds, or that lu,ls “re in [2*,+oo), and H— 
holds. Then there are one parameter families of upper and lower radial solutions 
with fast decay of denoted by frix) and iprix) respectively; hence tpriO) = 

D{t) = Cr(0) > 0, lim|,^|_>+oo |a;|"“^Cr(®) = Ac('^) lim|a,|^+oo \x\^~'^tjjrix) = 
L-ipir) > 0 where Lq{t) < The solution u{x,t;Q blow up in finite time, 

while the limit limt_>_|_oo ||u(a;, t; + |a;|)‘^||oo = 0 for any 0 < v < n — 2. 

Moreover ||Cr(a;)||oo = 1 |'*At(;c)||oo = D[t) —>• +oo, while Lq[t) < —>■ 0 as 

T — >■ — oo, while D[t) —>■ 0 and L.,i,{t) > L(^{t) +oo as t ^ +oo. 

1.7. Remark. For any fixed r G R, tfrix) > Ct{x) when x € M”. From the con¬ 
structive proof it follows also that both ||CT-(a;)(l-|-|x|‘^)||oo and ||'i/'T(a;)(l + |a;|‘^)||oo 
go to 0, as T —>■ -|-oo, for any 0 < v < m{ls), while they are uniformly positive for 
V = m{ls). 

A new aspect of Theorem 11.61 besides the generality of the potential we can 
deal with, is in the fact that we can find fast decaying initial data, with L'^-norm 
arbitrarily small, which blow up in finite time, while the critical decay indicated in 
literature (also by results as Theorem ll.5D for such a phenomenon seems to be slow 
decay, i.e. (see [13]). 

We emphasize that this result is new even when f{u, |a;|) = . Notice that, 

the dichotomy depicted in Theorem 11.61 and in Corollary 11.81 just below, takes 
place even for solutions slightly above or below a G.S. if we are in the hypotheses 
of Theorem fra The novelty here is that we can look at a much larger range of 
parameters and that this families of sub and super-solution have fast decay: Thus, 
we can find solutions with fast decay and L'^-norm small which blow up in finite 
time. 

The relevance of Theorem 11.61 follows from the next corollary. This latter result 
is an immediate consequence of the comparison principle. 

1.8. Corollary. Assume that we are under the hypotheses of Theorem \1.6\ Then 

for any e > 0 we can find smooth function (f : R" —> R, such that ||((>||oo < £ o,nd 
there is T^ > Q such that the classic solution u = u{x,t-,(j)) of U.l\) - U.2\} satisfies 
limt_>T’,^ ||u(a;, t; ^)||oo = +oo. On the other hand we can find smooth function 
(p : R" —>■ R, such that ||<(>||oo > l/e and the classic solution u = u{x,t-,<f>) of 11.11) - 
dll) IS dcfiTiGd foT QjTiy i ^ 0 CLTid sQ-tisfics lini ^—^_|_oq |x| )||oo — 0 

any 0 < u < n — 2. 
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From the above corollary we see how sensitive is, with respect to the initial 
data, equation We can find “large” initial data (j) which converge to 

the null solution and “small” initial data which blow up in finite time. Indeed, we 
can also construct initial data (pi and 4>2 such that for any e > 0 small we have 
||((^i — (^ 2 )(;c)[l + |a;|^]||oo < e, whenever 0 < v < m{ls), and u{x, t; pi) blows up in 
finite time, while u{x, t; P 2 ) is defined for any t and has the null solution as w-limit 
set. However we need to choose — p 2 ){x)[l + ||oo uniformly positive 

and bounded. 

Plan of the paper. In Section[5]we collect all the preliminary results concerning 
regular and singular solutions of m and, in particular, we prove new ordering 
properties. Section [3] is devoted to prove local existence of the solutions, in the 
classical, and in the mild case giving also a new result concerning singular solutions 
(which are slightly smaller than S.G.S. of (11.8^ 1. using a suitable weighted L°°- 
norm. Finally, in Section |3J we state and prove our main results on stability and 
long time behavior of the considered solutions. 


2. Ordering results and asymptotic estimates for the elliptic 

PROBLEM. 


The results of this sections, which are a key point for the whole argument, are 
obtained applying Fowler transformation to (11.81) . Thus, we set 


r = e 


( 2 . 1 ) 


m{l) = 


yi{s, 1) = 17(r)r’”« , y 2 (s, 1) = U' 

^ g(2/l,s;0 = /(2/le-™«^e*)e('"«+2)^ 


1-2 


Here and in the sequel I denotes a parameter which is always assumed to be larger 
than 2, so that m{l) > 0 (see the exemplifying case in (11.121) and also the parameters 
related to problem (12.31) below). Using this change of variables, we pass from (11.81) 
to the following system to which dynamical tools apply: 


( 2 . 2 ) 



m{l) 1 \ 

0 -[n-2-TO(/)] J \ y 2 


0 ^ 
-9{yi:S]l) J 


In the whole section the dot indicates differentiation with respect to s, and we 
introduce the following further notation which will be in force in this section: We 
write y{s, r; Q, I) = (?/i(s, r; Q; l),y 2 {s,T; Q)]l) for a trajectory of (12.21) where I = I, 
evaluated at s and departing from Q G at s = T. 

For illustrative purpose we assume first f{u,r) = so that we can set 

I = and system (Ea) reduces to the following autonomous system 


(2.3) 



m{l) 1 \ 

0 -[n-2-m{l)] J \ y2 


0 

-(2/1)^-' 


We stress that in this case we passed from a singular non-autonomous O.D.E. to an 
autonomous system from which the singularity has been removed. Also note that 
when (5 = 0 we can simply take I = q. 

System (12.3|) admits three critical points for Z > 2, = 2^(5^: The origin O = 
(0,0), P = {Pi,P 2 ) and —P, where P 2 = —m{l)Pi and Pi > 0. The origin is 
a saddle point and admits a one-dimensional stable manifold M and a one¬ 
dimensional unstable manifold M , see Figure [TJ The origin splits M (respec¬ 
tively M ) in two relatively open components: We denote by M® (resp. by M“) 
the component which leaves the origin and enters the semi-plane > 0. Since we 
are just interested on positive solutions we will call, with a little abuse of notation, 
M® and M“ unstable and stable manifold. 





To complete the depiction of the phase portrait in Figure [1] we recall the follow¬ 
ing result (see e.g. [7]) 


2.1. Remark. The critical point P of (ESI) is an unstable node for 2* < Z < cr,, an 
unstable focus if tr, < / < 2*, a center if ^ = 2*, a stable focus if 2* < ^ < cr* and a 
stable node if /> cr*, where 2», 2*, a* are as in EU and 


(2.4) 


n -2-1- 2^n — 1 

(Tjjc .- 2 ! — . 

n -f 2\/n — 1 — 4 


From some asymptotic estimate we deduce the following useful result (see, e.g. 
[1[5] for the proof in the p-Laplace context). 


2.2. Remark. Regular solutions u{r) of Equation correspond to trajectories 
Y{s;l) of system (12.2|) departing from points in M“ and viceversa. Positive solu¬ 
tions with fast decay u{r) of EH), correspond to trajectories Y(s; 1) of system (12.21) 
departing from points in M® and viceversa. 


Using the Pohozaev identity introduced in m and adapted to this context in 
[4] , we can draw a picture of the phase portrait of (12.31) , see Figure [U and deduce 
information on positive solutions of (11.81) : we postpone a sketch of the proof to the 
next subsection, where the general non-autonomous case is considered (anyway see 
[3] or [S] for a detailed proof in the more general p-Laplace context). Then it is 
easy to classify positive solutions: In the supercritical case {I > 2*) all the regular 
solutions are G.S. with slow decay, there is a unique S.G.S. with slow decay; in 
the critical case {I = 2*) all regular solutions are G.S. with fast decay and there 
are uncountably many S.G.S. with slow decay; in the subcritical case {2 < I < 2*) 
all the regular solutions are crossing, there are uncountably many S.G.S. with fast 
decay and a unique S.G.S. with slow decay. 

Since (ESD is autonomous we also get the following useful consequence. 

2.3. Remark. Fix U € M“ and S € M®. Gonsider the trajectories y{s,T-.,U), 
y{s,T\S) of (12.21) and the corresponding regular solution U{r,D) and fast decay 
solution V{r,L) of (11.81) . Then 

D = D{t) = £)(0)e-”"^ L = L{t) = L(0)e("-™)^ 

U{r, D) = DU{rD^/'^, 1) U(r, L) = Lv\rL^/^, 1) 

Proof. Since yi{s + = yi(s,0,Q) we get [/(re'^,Il(r))e"‘(*)'^ = U{r,D{0)) 

hence letting r —>■ 0 we find D{t) = ll(0)e“’"^*^'^, and this concludes the proof 
concerning U. Similarly we find 14(6®+"^, L(r))e"^®“'''^^e^™“"^'^ = U(e®, L(0))e"®, 
hence, letting s +oo we get T(r) = T(0)e("'“™^'^. Then again from j/i(s -b 
t,t,Q) = yi{s,0,Q) we get V{r,L) = LV{rL^^'^, 1): this concludes the proof. □ 


We stress that all the previous arguments concerning the autonomous Equation 
(EH) still hold true for any autonomous super-linear system (EH), more precisely 
whenever g{yi,s;l) = and g{yi\l) has the following property, denoted by 

GO (see [S] for a proof in the general p-Laplace context). We have the following 

GO: ■§^{yi', /) is a strictly increasing function for yi > 0 and 


lim 


Vi 


= 0 , 


lim 


g(i/i;0 

yi 


= -boo . 


In particular Remarks 12.1112.2112.31 continue to hold (see [S]). 

We emphasize that GO implies that is strictly increasing for yi > 0; then 

it follows easily that g{yi', 1) is strictly increasing too. 
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oh«< 1<2* 2*<l<aH« 

Figure 1. Sketches of the phase portrait of (EH), for g > 2 fixed. 


To draw correctly the analogous of Figure [T] for the present case, we need to use 
the Pohozaev identity introduced in [18] (see also [5] |3] for more details). Let us 
introduce the Pohozaev function 


P{u, u'; r) 




+ F{u,r)r", 


where F{u, r) = f{a,r)da. Now, consider the non-autonomous system (12.21) and 
denote by G{yi, s; 1) = g{a, s; l)da. In this dynamical setting the transposition 
of P(u, u'\r) is given by 

H{yi,y 2 ,s-J) = + G(?/i, s; 2*), 
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then, ify(s;2*) = (yi (5, 2*), 7 / 2 ( 5 , 2*)) solves (|2.2D with / = 2* we have the following 

J TT 

(2.5) —(2/1 (s, 2*), 2 / 2 ( 5 , 2*), s; 2*) = ^( 1 / 1 ( 5 , 2*), s; 2*) 

Moreover, if y(s; 2*) and 2 /(s; 0 trajectories of (12.211 corresponding to the same 
solution U{r) of (II. 8 L we get 

(2.6) H{y{s, 2*), s; 2*) = l),s; 1). 

where A{1) = n — 2 — 2m{l). We stress that (12.51) and (12.611 hold for the general 
non-autonomous system ( 1 ^ . 

Let us fix T G R and I > 2, and denote by 

'■= {(yi>2/2) I i?(2/i,2/2,T;0 = 6} 

K+ib) := {( 2 / 1 , 212 ) I H{yi,y 2 ,T;l) = b, 2/1 > 0} 


Then, there is 6*(t, /) < 0 such that the level sets K{b) of the function H, is empty 
for b < b*{T,l), they are two closed bounded curves contained in 2/2 < 0 < 2/1 and 
in 2/1 < 0 < 2/2 for 5* < 5 < 0 (the graph of the former gives K+{b)), K{b) is a 
8-shaped curve having the origin as center for 6 = 0, and it is a closed bounded 
curve surrounding the origin for 6 > 0. 

From ()2.5I1 we see that H{y{s, 2*), s; 2*) is increasing in s (respectively decreas¬ 
ing) along the trajectories y{s,2*) of (12.21) whenever G{yi,s;2*) is increasing in 
s (resp. decreasing in s). Moreover from (I2.6|l we see that H{y{s,2*),s;2*) 
and H{y{s,l),s;l) have the same sign. Thus, if we consider system (I2.3|l . for 
any Q G M“ and R G M® we get H{Q,s]l) < 0 < H{R,s]l) when I > 2*, 
H{R,s]l) < 0 < H{Q,s]l) when 2 < / < 2*, and H{Q,s;l) = 0 = H{R,s]l) 
ivhen I — 2 . Using (j2 .5 p and (1^ . it can be proved that the phase portrait of the 
autonomous system (12.31) is again depicted as in Fig. [U see e.g. mm- 

We collect here the values of several constants and parameters which will be 
relevant for the whole paper. Thus, recalling that m{l) = we introduce the 
followings 

(2.8) A{1) = n — 2 — 2m{l), C{1) = m{l)[n — 2 — m{l)]. 

Recall that {Pi,—m{l)Pi) is a critical point of (12.21) if it is s-independent, so Pi 
is the unique positive solution in y of gi{y,l) = C{l)y. When g{y,l) = y‘^~^ then 
Pi = (C'(O)^/^'?”^). Let n > 2 we denote by ct* < a* the real solutions of the 
equation in I given by 

(2.9) Ailf-4[C{l) + ^{Pi,l)]=0. 

which reduces to A{1)‘^ — 4((/ — 2)C{1) = 0 for g{y) = y'^~^. In this case the value 
of (T* coincide with the one given in (fm . 


2.1. The stationary problem: the spatial dependent case. Now we turn to 
consider (1^ in the s-dependent case. The first step is to extend invariant manifold 
theory to the non-autonomous setting; there are several ways to achieve the result: 
using skew-product semi-flow (see, e.g. [H]), or through Wazewski’s principle, see 
e.g. m- Here, we follow a simpler construction which is less general but preserve 
more properties (in particular the ordering results Propositions 12.8112. 9L used e.g. 
in [HIS]. So we introduce an extra variable, either z{s) = e^® or C(s) = e ™®, in 
order to deal with a 3-dimensional autonomous system. We use z and C in order 
to investigate the behavior respectively as s —?► — 00 (i.e. r —>■ 0), and as s —>■ +00 
(i.e. r -boo). 

We collect here below the assumptions used in the main results: 
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Gu: There is lu > 2* such that for any yi > 0 the function g{yi, s; lu) converges 
to a s-independent locally Lipschitz function g{yi, —oo; lu) ^ 0 as s ^ —oo, 
uniformly on compact intervals. The function g(yi,—oo]lu) satisfies GO. 
Moreover there is ro > 0 such that lims_)._oo 5 ( 2 / 1 , s;/„) = 0. Fur¬ 

thermore if lu =2*, we also assume that there is M > 0 such that g(yi, s; 2*) 
is monotone in s for for any j/i > 0 and any s < —M. 

Gs: There is Is > 2, such that for any 51 > 0 the function 5 ( 51 , s; Ig) converges 
to a s-independent locally Lipschitz function 5 ( 51 ; /g) ^ 0 as s —>■ -foo, uni¬ 
formly on compact intervals. The function g^°°{yi',ls) satisfies GO. More¬ 
over there is nj > 0 such that limg_>+oo -§^g{yi, s; Ig) = 0. Furthermore 
if Ig =2*, we also assume that there is M > 0 such that g{yi,s;2*) is mono¬ 
tone in s for for any ?/i > 0 and any s > M. 

A~: The function G(t/i,s;2*) := g(a,s;2*)da is decreasing in s for any 

2/1 > 0 strictly for some s. 

A~^: G{yi, s; 2*) is increasing in s for any 2/1 > 0 strictly for some s. 

Hypotheses Gu, Gs are used to construct unstable and stable manifolds for the 
Equation (12.21) when it depends on s, while A~ and A~^ mean that the system is 
respectively supercritical and subcritical with respect to 5 = 2 *, and are used to 
understand the position of these manifolds. 

2.4. Remark. Observe that if / is as in (HSl), dLll), (Ell) and (11.611 hold then Gu 
and Gs hold with lu and Ig defined as in (11.1311 and (11.1411 . 


Assume Gu. We introduce the following 3-dimensional autonomous system, 
obtained from (ESI) by adding the extra variable z = e™*: 

fyi\ (m{iu) 1 (y^\ ( ® \ 

(2.10) I 52 1 = ( 0 -[n-2-TO(;„)] 0 j ( 52 1 + ( - 5 ( 5 i,ii^;/„) 1 


Similarly if Gs is satisfied we set I = Ig and C{t) = e 

fyi\ (m{lg) 1 0 \ 

(2.11) h /2 = 0 -[n-2-m{lg)\ 0 52 

VC/ V 0 0 Vc, 


and we consider 


0 

+ I - 5 ( 2 / 1 ,- 
0 



The technical assumption at the end of Gu (and Gs) is needed in order to ensure 
that the system is smooth for z = 0 and C = 0 too. Consider (I2.10p (respectively 
(12.111) 1 each trajectory that may be continued for any s < 0 (resp. for any s > 0) 
is such that its a-limit set is contained in the z = 0 plane (resp. its w-limit set is 
contained in the C = 0 plane); moreover such a plane is invariant and the dynamics 
reduced to z = 0 (resp. C = 0) coincide with the one of the autonomous system 
(1121) where g{yi,s; lu) = 5 ( 51 , 5(yi, h) = 5 ( 51 , -boo; Ig)). 

Observe that the origin of (12.101) admits a 2-dimensional unstable manifold 
W“(2u) which is transversal to z = 0 (and a one dimensional stable manifold 
M® contained in 2 ; = 0), while the origin of (I2.11|l admits a 2-dimensional stable 
manifold W®(/s) which is transversal to the plane C = 0 (and a one dimensional 
unstable manifold contained in ^ = 0). Following [8], see also [151 0 we see 
that, for any r £ M, 

W“(t; lu) = W^ilu) n{z = , W®(t; Q = O {C = e—"} 

VF“(-oo; lu) = W^ilu) n {2 = 0} , W®(+oo; Ig) = W^lg) O {C = 0} 

are one-dimensional manifolds. Moreover they inherit the same smoothness as 
(12.101) and (12.111) . I.e., let AT be a segment which intersects 1 F“(to; lu) (respectively 
Is)) transversally in a point Q{tq) for tq £ [— 00 , -boo) (respectively for tq £ 
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(—oo, +oo]), then there is a neighborhood / of tq such that W^{t- ly) (respectively 
Vh®(r;/s)) intersects K in a, point Q{t) for any t G I, and Q{t) is as smooth as 
(12.1011 (resp. as (12.111) '). Since we need to compare and we 

introduce the manifolds: 

Is) := {R = gexp{ - - miQ^} gR'^\QG bb“(T; L)} 

W^iT-Jy) := {Q = i2exp{(m(/„) - m(l,))r} £ \R^w^iT-ls)} 

As in the s-independent case, we see that regular solutions correspond to trajectories 
in while fast decay solutions correspond to trajectories in Vh®, see [5J[S]. More 
precisely, from Lemma 3.5 in we get the following. 

2.5. Lemma. Consider the trajectory y{s,T,Q]lu) of i2.2\) with I = ly, the corre¬ 

sponding trajectory yit, r, R\ Ig) of ^2.2f) with I = Ig and let u(r) be the corresponding 
solution of ca). Then R = Qexp[{m{lg) — m{lu))T]. Assume Gu; then u(r) is a 
regular solution if and only if Q G or eguivalently R G W'^{t]Is). Assume 

Gs; then u{r) is a fast decay solution if and only if R G W^{T;lg) or equivalently 
QGW^T-ly)- 

Moreover if Q G bL®(r; ly) and ly > 2, then lims_>+oo y{s, r, Q; ly) = (0, 0), and if 
R G VL“(r; Ig) and Ig > 2 then lims_>_oo y(s, t, R; Ig) = (0,0). 

For the reader’s convenience we now report a result proved in which explains 
further the relationship between (11.81) and (12.21) : We recall that, close to the origin, 
bF“(T; lu) is locally a graph on the yi axis, while bF®(T; Ig) is locally a graph on its 
tangent space, i.e. the line y 2 = —{n — 2)yi. So let us consider a ball B{S) of radius 
(5 > 0 centered in the origin. Follow lF“(r;Z„) (respectively W‘^{t;Is)) from the 
origin towards j/i > 0. If d > 0 is small enough, we can choose a segment L C B{S), 
parallel to the 2/2 axis such that lF“(r;/„) (respectively VF®(t;Zs)) intersects L 
transversally a first time exactly in a point, say Q‘^(t) (resp. Q^ir)). We know 
that this point depends on r as smoothly as ( 1 ^ . so it is at least C^. 

Moreover, we have the following result analogous to 12.31 see BM- 

2 . 6 . Remark. Assume Gu. Consider y(s, t, (5“(t);/„) and the corresponding reg¬ 
ular solution U{r,a{T)) of (11.811 . Then a^r) —>• 0 as r —>■ —00 and a{T) —>• -l-oo as 

T —>■ -1-00. 

Similarly, assume Gs. Consider y(s, r, Q®'(t); and the corresponding fast 
decay solution V{r,l3{T)) of (11.81) . Then /3(r) —>■ 0 as t —>■ —00 and /3(r) —>■ -l-oo as 

T —>■ -1-00. 

Now we turn to consider singular and slow decay solutions of m- 
We observe that if ly > 2, then ( 12.1011 has a critical point in yi > 0, say {P 0), 
where P~°° = (P“°°, —m{ly)Pj~°°) is the critical point of the autonomous system 
(12.211 where g{yi,s-,ly) = g{yi,—oo]ly), and Pj~°° > 0. It is easy to check that 
(P“°°,0) admits an exponentially unstable manifold transversal to 2 ; = 0 which 
is 1 -dimensional (the graph of a trajectory which will be denoted by y“(s;/„)) if 
lu > 2*, and 3-dimensional if 2^. < ly < 2*. 

Analogously, if Ig > 2 * then ( 12.1111 has a critical point in j/i > 0, say (P+°°,0), 
where P+°° = (Pj^°°, —m{ls)Pi°°) is the critical point of the autonomous system 
( 12 . 21 ) where g{yi,s;ls) = g{yi,-\-oo]lg), and P^°° > 0. (P+°°,0) admits an expo¬ 
nentially stable manifold transversal to (^ = 0 which is 1 -dimensional (the graph of 
a trajectory which will be denoted by y®(s; Ig)) if 2^, < Ig < 2* and 3-dimensional if 
Ig >2*. In the whole paper we denote by f/(r, 00 ) the solution of (11.811 correspond¬ 
ing to y“(s;/„) and by y'^{s;lg) the corresponding trajectory of ( 12 . 21 ) with I = Ig-, 
similarly we denote by V{r,oo) the slow decay solution corresponding to y®(s;Zs) 
and by y®(s; ly) the corresponding trajectory of (I2.2|l with I = ly. 
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Note that if 2* < 4 < 2* < then the manifolds hF“(—oo; Z„) and M^®(+oo; Is) 
are paths connecting the origin respectively with p—°° and and contained 

in 2/2 < 0 < (we emphasize that this is not the case when 2^ < lu < 2* < Ig), see 
Figure [T] Using a connection argument we get the following. 

2.7. Remark. Assume Gu, Gs with 2, < /^ < 2* < then lU“(r; /„) and fU“(r; Is) 
are paths connecting the origin respectively with y'^{T;lu) and 2/“(r;Zs) for any 
r € [—oo,+oo); similarly W^{t;Iu) and W^{t;Is) are paths connecting the origin 
respectively with y®(r; Z„) and y®(T; U) for any r S (—oo, +oo] 

2.8. Remark. Assume Gu with lu >2*; then there is at least one singular solution 

U{r, oo) of (II. 8 |) . Moreover U (r, converges to Pi°° as r —0. Furthermore 

U(r,oo) is the unique singular solution if 

A specular argument gives us a similar condition for slow decay solutions. 

2.9. Remark. Assume Gs with Ig > 2*; then there is at least one slow decay solution 

V{r,oo) of (11.811 : moreover U(r, converges to P^°° as r —+oo. Such a 

solution is unique if 2 * < U < 2 *. 


Now we give a further result concerning separation properties which will be useful 
to construct sub and super-solutions for (11.81) . 


2.10. Remark. Assume Gs with Ig € [2*, a*) and consider two slow decay solutions 
U{r) and U{r) of (11.81) . Then U{r) — U{r) changes sign infinitely many times as 
r —> + 00 . Analogously, assume Gu with G (( 7 *, 2 *] and consider two singular 
solutions V{r) and V{r) of (11.81) : with a* and ct* in (11.91) and (12.41) . respectively. 
Then V{r) — V{r) changes sign indefinitely as r —^ 0. 


Proof. Denote by y(s) = y{s]lg), y{s) = y{s]lg) the solutions of (12.2|) correspond¬ 
ing to U(r) and U{r) respectively. Now assume Ig = 2* (and 5 ( 2 / 1 , s; 2*) monotone in 
s for s large). Then H{y{s), s; 2*) —>• b, and H{y{s), s; 2*) —>• 6 as s —>■ -boo, and 6 , b 
are both negative. lib>h>h* then y{s) converges to K+(b), and y{s) to K+{b), 
see (E3): by construction Kj^ ( 6 ) lies in the interior of the bounded set enclosed by 
Ar+( 6 ). Denote by A+ and A~ the point of K+(b) respectively with largest and 
smallest component 2 / 1 - When y{s) passes close to A+ we have yi(s) — yi{s) > 0, 
while when y{s) passes close to A~ we have 5 i(s) — 5 i(s) < 0, so the remark is 
proved. 

The argument works also if & > 6 = 6 *, so assume now b = b = b*, i.e. both y{s) 
and y{s) converge to 

We denote by h{s) = yi{s) — yi{s): note that h{s) —>■ 0 as s —>■ -boo and that it 
satisfies 


(2.13) 


h{s) + Bh{s) + N{h{s), s) = 0 


where B = B{lg) = ^{P+^)-milg)[n-2-m{lg)] = ^^(P+“)-iIi^ 


_ 


:>+oo\ {n-2) 


and 


N{h{s),s) :=g{yi{s) + h{s),s) - g{yi{s),s) --- {P^°°)h{s) = 


(2.14) 


=h{s) ( \^{yi{s) + (Th{s),s)- 
Jo ^oyi 


dyi 

dg+^ 

dyi 


.(P + OO) 


da 


So from ()2.14|) . Gs, and the fact that | 5 i(s)| -b |h(s)| —>■ as s —>■ -boo we see 

that N{h{s),s) = o{h{s)). Therefore for any e > 0 we find S = S{e) such that 
\N{h{s), s)| < e|h(s)| for any s > S. Note also that from Gs we get P > 0. Setting 


h{s) = p{s) 


cos{9{s)) 


(2.15) 


h{s) = p{s) sin( 0 (s)) 
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from (12.131) we get 


(2.16) 9{s) =—Vb — cos{9{s)) 


N{p{s) 


cos{9{s)) 

Vb ■■ 


s) 


pis) 


< -Vb{i 


e) < 


2 


for any s > S and S large enough. Since 9{s) —>■ —oo, and p{s) —>■ 0 as s —>■ +oo, 
but p{s) > 0 for any s S M, then h{s) changes sign indefinitely, and the Remark 
follows. 

Assume now Ig G (2*, tr*): then both y{s), y{s) converge exponentially to P+°°, 
therefore h{s) = y{s) — y{s) —0 as s —+oo. In this case (12.1311 is replaced by 


(2.17) h{s) - Ah{s) + Bh{s) + N{h{s),s) = 0 

where A = A{ls) = n — 2 — 2m{ls) < 0 and B = Bilg) = (P(*~°°) — m{ls)[n — 

2 — m{ls)] > 0. Note that \/B — ^ = 2e > 0 for /s G (cr,, 2*] and it equals to 0 for 
Is = <J*. So, using again (I2.15|) . and passing to polar coordinates we get 


(2.18) 


9{s) = -VB- 


Asm(26»(s)) 


— cos(0(s))- 


pis) 


A 


< -'/B + — +£<-£ 


So we find again that 9(s) —oo, and p(s) 0 as s +oo, thus h(s) changes 
sign indefinitely, and the Remark follows. 

The case of singular solutions V(r) and V (r) can be obtained from the previous 
repeating the argument but reversing the direction of s. □ 


Following we can show that if A~ holds then (11.81) is supercritical, while if 
A+ holds then (11.81) is subcritical. To be more precise we have the following (see 
0 Theorems 4.2 and 4.3]). 

2.11. Proposition. Assume Gu, Gs, with ls,lu G (2*, 2*], and A+, then all 
the regular solutions U(r,a) are crossing, i.e. there is R{a) such that U{r,a) > 0 
for 0 < r < R{a) and U{R{a),a) = 0. Furthermore R{a) is continuous and 
R{a) +c» as a —>■ 0, and if lu < 2* then R{a) —?► 0 as a ^ +c». 

Moreover, all the fast and slow decay solutions are S.G.S. So for any /3 > 0 the 
fast decay solution V(r,j3) is a S.G.S. with fast decay; if Ig < 2* there is a unigue 
S.G.S. with slow decay, say V[r,oo), while if Ig = 2* there are uncountably many 
S.G.S. with slow decay. 

Proof. This result is borrowed from O Theorem 4.2] , where it is proved in the p- 
Laplace context in a more general framework, so here we just sketch the proof. The 
main idea is to use the Pohozaev identity as done in the previous subsection: From 
(12.51) we know that the function H{y, s; 2*) is decreasing along the trajectories, and 
it is null for y = 0. Using also (12.61) we see that if Q G IU“(t; 1^) and R G IU®(r; Ig) 
we get H{Q,t-,Iu) > 0 > H{R,T;lg). Recalling which is the form of the level 
set K{b) of H (see (12.71) and the discussion just after it) we deduce which is the 
position of kF“(r; lu) and kF'’(r; Ig) and using Lemma E31 Remark [^?71 we conclude 
the proof. □ 


With a specular argument we get the following. 

2.12. Proposition. Assume Gu, Gs with lu,lg > 2*, and A~, then all the regular 
solutions U{r,a) are G.S. with slow decay. Moreover all the fast decay solutions 
V(r,j3) have a positive non-degenerate zero r = i?(/3), i.e. V(r,j3) is positive for 
any r > i?(/3) and it is null for r = R{(5). Furthermore R{P) is continuous and 
R{P) 0 as ft ^ +CXD, and if lu > 2* then R{ft) +c» as ft ^ Q. Further if 
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Z„ = 2* there are uncountably many S.G.S. with slow decay, while if lu > 2* then 
there is a unique S.G.S. with slow decay say U(r,oo). 

Now we give a Lemma, consequence of Propositions 12 .1 ll and 12.1^ which allows 
to extend picture [T] to the non-autonomous setting. Assume A~h,Gu,Gs with 
2 < lu < 2* and 2, < 4 < 2*. Follow W'^{t;Iu) from the origin towards := 
{(yi, 1 / 2 ) I yi > 0 }: it intersects the t /2 positive semi-axis in a point, say We 

denote by the branch of W^{t-,Iu) between the origin and and 

by E^{t) the bounded set enclosed by W^{t; lu) and the 2/2 axis. 

Similarly assume A~,Gu,Gs with lu > 2* and Ig > 2*. Follow W^{t]Is) from 
the origin towards yi > 0: it intersects the 2/2 negative semi-axis in a point, say 
(3®(r). We denote by W^{t-,Is) the branch of between the origin and 

Q®(r), and by E^{t) the bounded set enclosed by lF“(r; lu) and the 2/2 axis. Using 
the fact that H{Q,t;Iu) > 0 > H{R,t;Iu) for any Q S LF“(r;/„), i? G LF®(r;Z„) 
if A+ holds and 2 < /„ < 2* and 2* < < 2*, while E[{Q,t;Is) < 0 < H{R,t;Is) 

for any Q G LF“(r;/s), R S W^{t;Is) if A~ holds and lu > 2 *, Ig > 2 * we get 
the following Lemma, which is useful to construct a new family of sub and super¬ 
solutions, see also Remark [521 

2.13. Lemma. Assume A~h, Gu, Gs with 2* < /„ < 2* and 2* < /s < 2*. Then for 
any r G K, W‘^{t;Iu) C E'^{t); assume further Ig < 2*, then W’^{t\Iu) is a path 
joining the origin and 2 /®(t;Z„). 

Assume A~, Gu, Gs with lu > 2* and Ig >2*. Then for any r G K, LF“(r; Ig) C 
E’^{t); assume further lu > 2*, then LF“(t;Zs) is a path joining the origin and 
y^{T;lg). 

We emphasize that the sets E‘^{t), E^{t) have the following property: let Q G 
E'^{t), R G E‘^{t), then y{s,T,Q-,lu) G E'^ft) for any s < t, and y{s,T, R-,ls) G 
for any s > r. 

When lu = Ig = 2* we have a slightly different situation. Denote by P*{t) = 
{P* (t), Pf (t)) the critical point of the autonomous system (12.21) where I = 2* and 
22 ( 2 / 1 , 3 ; 2*) = 22 ( 2 / 1 , t; 2*). Denote by Pf := inf{Pf(T) | r G R}; in this setting we 
have Pf > 0 and we denote by P* = Pf /2. We denote by P* = (P*, —m(2*)P*). 

2.14. Lemma. Assume Gu,Gs with lu = Ig = 2*. Assume further A^, then 

for any r G M the line 2/1 = P* intersect the manifold IF“(t) in = 

(P*, (^ 2 and in Q'^’~{t) = (P*, ('’"))> and it intersects W‘^{t) in {t) = 

{P*,Q 2 ~{t)) and Q 2 ’~ (t) < Q 2 ~ < —m{2*)P* < Q^’^ir). Moreover, if y{s) 

corresponds to a S.G.S. with slow decay, there is Q = (Qi,Q 2 ) S lF“(r) such that 
Qi = yi{T) and Q 2 > 2 / 2 (t) for any r G K. 

Now, assume A~, then for any r G R. the line 2/1 = P* intersect the manifold 
W%t) in Q®’±(t) = {P*,Ql'^iT)), and VF“(r) in Q^’+{t) = (P*,Q“’+(r)) and 
Q 2 ~ (t) < Q 2 ’~ i^) < —m(2*)P* < Q 2 '^{t). Moreover if y{s) corresponds to a 
S.G.S. with slow decay, there is Q = {Qi,Q 2 ) G W’^ir) such that Qi = 2/1 (t) and 
Q 2 < 2/2 (t) for any r G R. 

Proof. We recall that LF“(t; 2*) and LF'’(r;2*) depend smoothly on r and that 
they become the graph of a homoclinic trajectory as r —> —00 and as r —>• -|-oo 
respectively. Denote by 

5'(t) := {(yi, 2 / 2 ) I H{yi,y 2 ,T; 2 *) = 0 , 2/1 > 0 } 

and by H'^{t) = (P*,iL+(r)), and H~ (t) = {P*, T[~ (t)) the intersection of S{t) 
with the line 2/1 = P*, where T[~ (r) < P+(r). 

From an analysis of the phase portrait relying on Wazewski’s principle it fol¬ 
lows that 1 F“(t; 2 *) (respectively 1 F'*(t; 2 *)) intersects the line 2/1 = P* for any 
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Figure 2. Sketch of the proof of Lemma [2.131 when A+ holds. 


T G R, see [B] for a proof in the p-Laplace context. Follow VF“(r; 2*) and 1 F®(t; 2*) 
from the origin towards yi > 0: we denote by the first intersection of 

bF“(T;2*) (resp. of bF®(r;2*)) with the line yi = P*, and by the first 

intersection of 2*) with yi = P*. Using transversal smoothness of the man¬ 

ifold lF“(r;2*) and VF®(t; 2*), see subsection 2.1, we see that we have at least a 
further intersection with such a line, respectively for r << 0 and for r >> 0. We 
denote by the second intersection of VF“(t„; 2*) with the line yi = P* 

and by the second intersection of VF®(ts; 2*) with the line yi = P*, for 

any < —N and Ts > N and N > 0 large enough. Set = {P*, {tu)) 

and (5®’^(rs) = {P* , (ts))'- possibly choosing a larger N we can assume w.l.o.g. 

that Q 2 ’^{tu) > -m{2*)P* > Q 2 ~{tu) and Q 2 ^{tu) > -m{2*)P* > 

We denote by VF“(tu) the branch of 1F“(t„;2*) between the origin and 
and by W‘‘{ts) the branch of W^{ts] 2*) between the origin and Q®’+(rs). 

Assume A+; then lF“(r; 2*) lies in the exterior of the bounded set enclosed by 
S{t) for any r: we claim that exists for any r G R. In fact consider the 

semi-line L{t) = {(P*,i/ 2 ) | J /2 < the flow of (12.211 on L{t) points towards 

?/i < 0 for any r G R. Hence the trajectory y{s,—N,Q'^’~{—N)]2*) crosses the 
line yi = P* for s = —N and then the yi < 0 semi-plane, and similarly for any 
Q G lU“(r) the trajectory y{s, —N, Q, 2*) will cross the line yi = P* for a certain 
s > —N and then the yi < 0 semi-plane. Hence, for any r > —N, the branch of the 
manifold IF“(r; 2*) between the origin and y(r, —N, Q'^’~{—N); 2*) will surround 
S (r) untill it crosses a second time the line yi = P* and the claim is proved, so we 
get pictured 

Now denote by D^{t) the bounded set enclosed by IF“(r), the segment between 
and H~{t) and the branch of S'(r) between H~{t) and the origin; ob¬ 
serve that by construction if Q G D^{t), then y{s, r, Q, 2*) G D^{s) for any s < t. 
Since S'(r) C D‘^{t) we see that if y{s) corresponds to a S.G.S. with slow decay, 
then y(s) G S{s) for any s G R. 

Reasoning in the same way but reversing the direction of s we see that if A~ holds 
then we can construct W^{t) for any r G R. Denote by D‘’{t) the bounded set en¬ 
closed by the segment between and and the branch of S{t) 

between H+{t) and the origin. Then if y{s) corresponds to a singular solution, 
then y{s) G S{s) for any s G R. So Lemma [2.131 follows. □ 

Now we give a Lemma which is useful to detect the w-limit set of solutions of 
( 11 . 11 ) ( 11 . 21 ) in the case where </> is a radial upper or lower solution of (11.81) . 
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2.15. Lemma. Let U{r) and V{r) be positive solutions of il.8\) either regular or 
singular and assume that there is Z > 0 sueh that U(Z) = V{Z) and U'{Z) < 
V'{Z). Denote by 


(2.19) 


r V{r) ifr<Z ^ < Z 

\ U{r) ifr>Z ’ ’ \ V{r) if r > Z 


Assume A~,Gu,Gs with > 2* and Is £ [2*,ct*). 

Then U.8\) admits no solutions (j)(r) either regular or singular such that 0 < 
^ C(^) solutions (j)(r) such that <f)(r) > ipir) for any r > 0 . 


Proof. From Proposition 12.121 we know that all the positive solutions have slow 
decay. Assume first Is £ [ 2 *,( 7 *). Then from Remark 12.101 all the slow decay 
solutions of cross each other indefinitely as r —> +oo, so the Lemma easily 
follows. □ 


Reasoning in the same way we get the following: 

2.16. Lemma. Let U{r), V{r), ((r) and tpir) be as in Lemma \2.15\ Assume 
A.+ , Gu, Gs with 1^ £ ((T», 2*] and Is £ (2,; 2*]. 

Then U.8\) admits no solutions (j)(r) either regular or singular such that 0 < 
^ C(^) solutions (j){r) such that (f(r) > ip{r) for any r > 0. 


3. Local existence 


In this section we introduce some basic facts and definitions related to the prob¬ 
lem (ll.ll) - (ll. 2 l) . and exploiting techniques similar to those used in [571 §lj § 2 ] (see 
also Ch. II]), we prove local existence for the solutions of problem (fmi-lfTTI). 
For the remainder of this section we will make the following assumptions, in addi¬ 
tion to FO, on the potential f in (HH), i.e.: 

Fu: Gu holds and there is H > 0 such that -^[g{yi, s]lu)\ < D\yi\^ for any 
0 < 1/1 < 1 , and 0 < s < 1 . 

Fs: There are ^>0, C'>0, (5>0, e>0 such that \f{u,r) — f{u h,r)\ < 
Chu^r^ whenever r > 1 , 0 < u < 1 , and 0 < h < e. 

Assumption Fu is very close to Gu (and it is actually satisfied in all the motivating 
examples given in the introduction), while Fs is more standard and it is adapted 
from m- Let us introduce the following map 

(3.1) w(x) = I j^j ^/,5 where 0 < 1 /< m(Z„). 

We emphasize that ii i = 0 then w{x) = 1 for |a;| > 1. Moreover if we set i/ = 0 
then w{x) = 1 for |a:| < 1 so we are dealing with bounded solutions, while if we set 
n > 0 we can deal with solutions which are unbounded for \x\ small and are not 
defined for x = 0 . 

Let us recall the definitions of continuous weak solution and Cs-mild solution 
to the problem (HID-dia). 


3.1. Definition. We say that a function m is a continuous weak (c.w.) solution of 
(ini) 1121) if u is continuous and it is a distributional solution: i.e. if u{x, 0 ) = (j){x) 
and, for any y £ C^’^(]R" x [0,T]) with y > 0 and supp? 7 (-,t) (s K" for all t £ [0,r], 
it holds true that 


(3.2) 


u{x, s)ri{x, s)dx\Q^ = 


Ti 


u{x,s){yt + Ay){x,s) + f{u, |a;|) 77 (x,s) 


dxds 


jk" Jo Jr^'- 

A Ti £ [0,r]. Further, u is a c.w. lower (respectively upper) solution of (ll.ll) - (ll.2D 
if u(x,0) > (f(x) (resp. u(x,0) < 4>{x)) and we replace “ = ” in (13.211 by “ > ” 
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(resp. by “ < We call a function u a classical solution if it satisfies (ll.ll) - (ll.2l) 
andue x (o,r))nc(R" X [0,r]). 


'■ oo /TtLn 


(3.3) 


Let (j) G := n L 

:= (47rt)“"/^ f exp ( 

Jr”- ^ 

i-t 


). We introduce the following operators 
X - yp' 


4t 


-)Hy)dy- 


F^{u) = + J ed I • |)ds)(x) 


3.2. Definition. We say that u is a Cs-mild solution of (ll.l|) - (ll.2|) on R" x [0,T) if 

(3.4) u G Cb(M’" X [0,T']) := C'(]R” x [0,T']) n x [0,r']), 0 < T' < T, 

(3.5) u{x,t) = F^{u{xfi)), for(x,t) G R" x [0,r). 


Also, we say that u is a Cs-mild lower solution (or upper) solution if “ = ” in (13.511 
is replaced by “ > ” (“ < ” respectively). 


We also consider the norm ||^||x := ||'('w^||l°°(R") and the weighted space X = 
L“(R”) = {(j) I ||(('||jc < +oo} where w is defined in (13.11) . We denote by Cs(R" x 
[0, T']) := C((R"\ {0}) X [0, T'])nL“([0, T'], A) and we give the following definiton 

3.3. Definition. We say that u is a Cg-mild solution to ((ri|)-(lL2l) on (R"\ {0}) X 
[0,T) if u G (^^(R" X [0,T']) for any 0 < T' < T and satisfies u{x,t) = F^{u{x,t)) 
for (x,t) G (R"\ {0}) X [0,r). 


Note that if > 0 and u is a Cs-mild solution then it may be unbounded as 
a; —^ 0. Therefore we can deal with initial data and solution having a singularity in 
the origin, and we will prove local existence and uniqueness for such initial data. It 
is worth recalling that with our assumptions we have singular stationary solutions 
fis which behave like as a; —>• 0. 

An interesting question, still open even for the starting case /(u) = is 

whether stationary S.G.S. are stable or not. One of the difficulties is in fact to 
prove local existence and uniqueness for nearby initial data (which is in general 
violated but maybe recovered in some special space and for some parameters, see 
[25]!. In fact we cannot even hope for a general local uniqueness result: We need to 
prescribe a class of function within local uniqueness is recovered, due to the presence 
of self-similar solutions converging to singular data (see, e.g. [551 US), and to a new 
class of solutions with moving singularity recently described in [23l [24] . 

We stress that here we are forced to stay below fis since we need to require 
v < m{lu), so we cannot start a stability analysis for stationary singular solutions. 

The following result is a direct consequence of m Lemma 1.5]. 

3.4. Lemma. Assume Fu, Fs. Let u be a continuous weak upper (lower) solution 

of with n > 3. Assume that there exist k,(3 > 0, and 0 < a < 2 

such that f{u,x) < fcexp(/3|a:|“) on R" x [0,T]. Then u{x,t) > {<)F^{u{x,t)) on 
R” X [0,T]. 

3.5. Remark. By this lemma it follows that a c.w. solution of (HID-dra) satisfying 
either (13.41) in Definition 13.21 or the analogous weighted condition in Definition 13.31 
is also respectively either a Cs-mild solution or a Cs solution. The converse is also 
true (see the proof of [271 Lemmal.5]). 

To prove Lemma [331 it is sufficient to adapt the proof of [571 Lemma 1.5] to the 
present case. By Lemma 13.41 we also have the next result. 


3.6. Proposition. Assume Fu, Fs. Suppose that u is a continuous weak upper 
(resp. lower) solution of in R"\{0} such that ||m||x is bounded, then u is 
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a Cs-mild upper (resp. lower) solution of (11.11) - (11.21) . The converse is also true, 
provided 4){x) > u(x, 0) (resp. (j){x) < u{x, 0) ). In particular if we set v = 0 we see 
that a continuous bounded weak upper (resp. lower) solution of o is a Cs-mild 
solution of (fLT])-(fT:2]) and viceversa 

Take p > 0 and f) G X, and denote by Bp := Bj (0) C L°“([0, T\,X) the ball of 
center 0 and radius p, in X; the radius p will be chosen properly later. We now 
prove local existence and uniqueness for Cs and Cg-mild solutions. 

3.7. Lemma. Assume Fu, Fs. If the initial datum (f G X, there is > 0 such 
that the operator F^{u) defined by FI A) has a unique fixed point in i?J (f)) for any 
0 < T < Tff,, and ifT,j, < +oo, then lim^ ||u(-,t)||x = +oo. 

<f> 

Proof. We claim that the operator maps Bp in itself and it is a contraction: 
Then the Banach fixed point theorem provides existence and uniqueness of a fixed 
point u for F^. 

Observe first that if u G Bp, then \u{x, t-,(j))\ < pw{x)~^ a.e. in M" x [0, T]. Here, 
we take p = 2{2Di + 2*"+^ + 2^/‘^)||(/)||jf, where Di := 

From relation Fu, for \x\ < 1, we get 

f[u, |a:|) < f{pu:-\x), IccI) < 5(p|xr('“)-Mn(|x|); 


where k = Dp^^^. Then, for any v G Bp, we also have 

r-i 9f 


\f{u,\x\) - f{v,\x\)\ = 

dg 


f df 

J — (sM + (1 - s)u), |x|)[m - ujds 


(3.7) 




where we have redefined k~ = Hmax{p^+^, p'^}. From Fs, for |a:| > 1, we also get 
(3.8) f{u, |a;|) = f{u, |a:|) - /(O, |a:|) < C\u\^+^\x\^ < Cp^+^\x\~^/^ < k+w{x)~'^, 

(3 9 ) < C-^\\u-v\\^\x\^ 

< k~^w{x)~^\\u - u||x, 

with fc+ =C'max{p^+^, p^}. Till the end of the proof we need the following straight¬ 
forward estimate: Let A G (0, n) and denote by F the Euler Gamma function, then: 


(3.10) 


e 4 


-dp < 


/R" 


(47r)"''^|77|^ 2 " i(n-H)r(n/ 2 ) 


+ e ^ dp < C{A). 


and we can set C{A) = 2 if H G (0, n — 1). We now proceed to prove that F^f,: Bp 
Bp. From (13.31) we have that 


\F<l>{u)\{x,t) < \e*^(j)\{x) 


“■>(-|3S) 


0 Js." (47r(t — s))”^^ 
< w~^{x)\\e*^(j}\\x(t) -\-w~^{x)I, 


fiu{y,s),\y\)dyds 


(3.11) 
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where 


I = w{x) 


exp - 




/O \J\y\<l •l\v\>l. 

Using (I3.6|) - (I3.7I) we get the following 


(47rs) 


^/2 


/(■w, \y\)dyds =: 4 + h- 


(3.12) 


la <k 


■/ 

-'ll/I 


exp ( - 


0 V-/lfl<|y|<l •^\y\<^ ^ u)+i' 

=:k-{I-+I+) 


dyds 


and 


I- < 2’' 


exp - 






0 -/lfl<|y|<l (47rs)”^^|?/P S(m{la) v) 


dyds 


0 Jwx- (47rs^ ^ |yP 

exp(-|77p)|77|'5W'“)-U- 


^^77 

Jo Jr" 


(7r)ra/2gl-5(m(Z„)-i/)/2 


dyds 


'-dyds < 


where iiT” > 0 is a positive constant, and we used the fact that the convolution of 
radial decreasing function is radial decreasing too (see [571 Lemma 1.4]), and that 
n — 3 + S{m{lu) — 1 ^) > —1. 

Since ja; — yj > ||a;| — |y|| we get 


/+ < 


< 


< 


w{x) exp ( - 


-dyds 


0 (47rs)"^^|2/|2+‘^ 

w{x) exp ( - exp ( - i|f) 


0 Ay\<^ 1^) 


dyds 


Jo 


w{x) exp ( - .yj) 


-[S{m{l^)-v)/2\+v/2 


ds 


I 


exp - 


(-¥) 


< 2”/2C' 


7(x) 


(47r)"' |7^|2+i^-<5(m(iu)-!/) 


dy 


gl — Sm{lu ) l2-\-y j2 


exp 


V 32s/ 


ds 


where we used that 2v — 5{m{lu) — v) <2-\- m{ly) < n, and C = C{2 + m(l„)) 
is the constant defined in (13.1011 . Now we need to distinguish between the |a;| < 1 
and the |a;| > 1 case. Assume the former so that w(x) = |a;|‘", and observe that 
h{a) := < Di = h{\Qv)-, for any a > 0 we get 


(3 13) /+ < ZliC f -- ds < 

^ ^ 1 '-' gl-[5(m(i„)-i/)/2]“* - 

where = ^(m(7)-t/)/2 > 0 is a constant. When |a;| > 1 so that w{x) = 
setting h{a) := we find h{a) < D 2 := h{32i/6) for any a > 0 and 
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similarly h{a) := e-“/64ai-[-5(™(i.)-i'-i/^)/2]+^/2 jg bounded, say h{a) < D 3 . Thus 


I+<C 


(3.14) 


/o 


e 

2 \ 2d 


exp - 


( 

V 64s/ 


exp 


(-») 


gl—[6{m{lu) — i^— 1 /^) l2]-\-y 12 


ds 


ft _ |a;|2 _ 1 

< CD 2 / h('-^)h(-)ds < CD 2 D 3 t < K+t 
Jq s s 


with a possibly larger constant . Now we estimate h- From (IHli-dlll) we get: 

^ _ \x-v\^ 


(3.15) 


h < k~^w{x) 


exp I 


4s 


dyds 


lo V /|y|>max{^,l} JJ (47rs)"'^^|?/|^/^ 
=:A:+(4-+/+). 

Observing that h{a) = a^/(2-5)e-“/32 < hi^U/d) and Ch = 2^'^/^h{8i/S), we find 


h <w 


Jo J\v\ 


exp - 


(-4#) 


0 /|i/|>max{Jfi,i} (47rs)"^^|y|^/'5 " Jo w(max{f|f, 1}) 


dyds < 


j{x) 


kl 




\x\^ 

It < I w{x)e 32s ds 


_i£r 

e 32s 


jdy < 2^"/2|a;|^/^te"^ < , 


Ji<\y\<'4 (47rs)”^^ 

Therefore, there is A" > 0 such that / < ATmax{t^+^/*^^^\t,with 
AT = K{n,£,iy, ||<(>||x) and relation (13.1111 reduces to 

(3.16) \F^{u)\{x,t) < w-i(x)||e‘^0||x+w-i(a;)A:max{ri+^/(2O^r^7’5(™(A)-O/2}_ 

To estimate the term (x)||e^*(/)||x we follow an approach similar to the one 
used above. Indeed, we rewrite e*^(j){x) as follows 


(3.17) e*^(l){x) = 

Hence, we get 

la < Wt^Wx 


'|y|<i /|y|> 


exp 


(-1^) 


(47rt)’ 


-j^-^(t){y)dy =: /« + Ip. 


exp 


(-‘^) 


dy =■ da +It- 


For ^ < bl < 1, we reach 


IF < 


exp - 


(-4#) 


(3.18) 


\mx2'^ f 

w{x) Jl^<\y\<i 


dy 


“ w{x) 


exp 


/R" 


(ll) 

n/2 


(47rt) 


dy = 


Mx2^ 

w{x) 


For the term /+, using (13.1011 we obtain 


It < ll</'IUexp(- 


exP ( - ¥) 

32f / 7|y|<min{lfi,l} (47rt)"^^|?/|^ 


dy 
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Now, arguing as in (I3.13L (13.141) we find for |x| < 1 

1//2 


(3.19) 


It < m\x 


(¥)' exp(^) ^ 274111^1 




la;!*" w{x) 

while for |x| > 1, since h{a) = a^/(2'5)g-a/64 ^ gg|. 

X/(2S) 


it<m\x 


exp exp 2 D 2 M 


w[x] 


iy!2-ll{25) 


< e(t) 


X 


w{x) 


where £{t) := exp —>■ 0 as f —>■ 0. So choosing t small enough we 

can assume that (13.191) holds for |a;| > 1, too. Take into account 7^, to get 


h < 


ll'/'ll 




(3.20) 


_ 

exp ( - ) 


^ ll^llx II^11 


(4ir) 


i/2 


dr/ 


< < &max{2^2^/n 

w{x/2) w{x) 

Collecting the estimates (I3.17I) - (IXT^ - (I3.191) - (I3.20I) . we have that relation (I3.16p . 
and hence (13.111) . for T < 1 gives 

||F,^(u)|U(T) < (277i + 2"+i + + 77max{r, 

< p/2 + C'T||()>||x +77max{T,T^('"('“)-‘')/2|^ 

Let To = To{n,£, ||(/)||x, p) > 0 be such that, for any T < Tg, 

77max{r,r‘5(™('“)-'^)/2} < ^/2. 

Then, we have that 

II^<^{u)||l=o(o,t;X) < P, for any T < Tg 
and hence maps Bp into Bp, for T <Tq. 

Analogously, let u,v £ Bp, we get 

fc- exp - -i^^) 


(3.21) 


\F^{u) - F^{v)\{x,t) < f (f 
Jo ^J\i 


0 '"e/|y|<l ^47rs)"'^^|y|2+‘2-5(m(Z„)-i/) 

k+ exp - -i^^) 


dy+ 


dyj\\u- v\\xds. 


Repeating the argument of (13.121) and (13.131) we get 

(3.22) ||[F,^(u)-F,^(u)]|U(t) < fc-[i7-t^(’"('fo-fo/2 + 7C+t] + 7+[2^/^ + at^/(25)]^^ 


Therefore, taking T < Tq sufficiently small, it follows that F^ maps Bp into Bp and 
it is actually a contraction. From the contraction principle, we obtain existence and 
uniqueness of a fixed point u in Bp which in turn implies the existence and local 
uniqueness of a C^-mild solution to dnii-doi). Then, we can restart the reasoning, 
by setting (p{x) = u{x,T) and go up to T,p by a ladder argument. Note that if 
< 00 and limt-pr^ ||u(x,t)||x is bounded we can restart the ladder argument 
and obtain a continuation interval [0,T'] D [0,T,^) and this is a contradiction. 
Hence if < c» we get limt_pT,^ ||it(a:, t)||j(: = + 00 . □ 
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As a direct consequence of Lemma 13.71 we have the following existence result 

3.8. Theorem. Assume Fu, Fs. Let (j) G X he the initial datum for the Equation 

(HH). Then problem inii-dia) has a unique weak solution u on R” x [0, T^), > 

0. Iff) G Cb{M.") = C'(R")nL°°(R”) thenuG C'(R’" x [O,T0)) x [0,T^)), 

and if < +oo then limj_j^j,- ||u(-, t)||ioo(Rn) = +oo. Similarly, if f) G C's(R") = 

C'(R"\{0})nL“(R”) t/ien M G C'((R"\{0}) X [0, T,^)) n L'^([0, T,^), L“(R’^)), and 
ifT^p < + 00 , thenlim^^j.- ||■u(•,t)u;(-)||L-=(R") = +oo. 

Furthermore, */</>> 0, then u > 0; if f) is radial, then u is radial in x; if f) is 
radial and radially non-increasing, then u is non-increasing in |x|. 

3.9. Remark. Assume that / is locally Holder continuous in x and assume also there 
exists I > 0 such that f{u, |a;|)|a;|“* = f{u, r)r~‘ is locally Lipschitz in u uniformly 
with respect to x (and so in r) in any bounded subset of R". In such a case, using 
[271 Lemma 1.2] and arguing as in [27], one can verify that the Cs-solutions are 
actually classical. 

In particular, we have that the potentials in (IT2T1)-(fT31) verify these conditions. 
As a consequence of Proposition lS.bl and Lemma 13.71 we have the following result 

3.10. Theorem. Assume that Fu, Fs are verified. Then 

(i): Suppose that u and u are Cs-mild upper and lower solutions of (HID-lIll) 
on M" X [0, T). Then u > u on R” x [0, T), and the unique Cs-mild solution 
of (ll.ll) - (ll.2l) on R" X [0,r,^) satisfies that u < u < u on R" x [0,T) and 
70 > T. 

(ii): If the initial value f) in (na) is a c.w. upper (lower) solution of ra, then 
the Cs-mild solution u of (frT])-(fT:2]) is non-increasing (non-decreasing) in 
t G [0, T^). 

(Hi): If (j) is radial, then u, u and u are radial for any t in their dominion of 
definition. 

(iv): If cj) is a c.w. upper (lower) solution but not a solution of |7.7| ), then 
ut{x, t) < (>)0, t > 0. 

The proof of this theorem is omitted because it can be easily derived by adapting 
that one of [271 Theorem 2.4] to the current case. We point out that claim (iv) 
follows directly by exploiting a comparison principle, arguing as in [271 Lemma 2.6] 
(see, e.g., m for a full-fledged proof of this well-established comparison argument. 
See also [2]). Further, a result analogous to Theorem 13.101 holds true also in the 
even simpler case of the Cs-mild solutions of (dllll-dllll) on R” X [0,r) (see [27] 
Theorem 2.4]). 

To conclude this section we give a result about the global solution u{x,t;(j)) of 
the problem (ll.lll - (ll.2L for f) G Cs(R") or f) G C 5 (R"). 

3.11. Remark. Assume that (pis a singular upper (respectively lower) solution. From 
Theorem 13.101 point (iv), which translates (27] Lemma 2.6] to the current case, it 
follows that limt_>T,^ u{x, t; cp) = u{x, T^; p) exists for any x ^ 0. Following the proof 
of Claim 2 of (27] Theorem 3.6], using Lebesgue dominated convergence theorem 
and regularity theory for elliptic equation we see that u{x, T^j,] p) is a distributional 
solution of (da. Moreover if p is radial then u(x,T^\p) is radial too. 

4. Long time behavior: main results 

Now we are ready to state and prove our results in their general form, from 
which Theorems and Corollarv ll. 81 follow. 

Let w{x) be defined as in (EH). 
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4.1. Theorem. Assume Fu, Fs, A , Gu and Gs, with 2* < 4 < o'* lu > 2*. 

(i) : If 4>{x) ^ U{\x\,a) for some a > 0, then it holds that ||it(x, (^)||oo —>■ 0 

and ||m(x, t; (/))(! + |a;|‘^)||oo —>-0 as t ^ +oo for any 0 < i/ < m(ls). 

(ii) : Let (p be a continuous initial data or a singular one such that ||(/)(x)4(;(x)||oo 

is bounded, for some v G [0,m(/s)). If <f)(x) ^ U{\x\,a) for some a > 0, 
then \\u{x,t] (j))w{x)\\oo must blow up infinite time. 

This result establishes that G.S. -with s.d. are on the border between the basin 
of attraction of the null solutions and initial data which blow up in finite time, in 
the range of the considered parameters. We emphasize that we need Ig < <J*, but 
lu has no upper bound: this allows e.g. f{u, r) = (1 + r°')u‘^ where q > a* and 
a G {—^{q — O'*), ^^(9 — 2*)). It is easy to check that A~ might be replaced by 
H—, and that Theorem 14.11 implies Theorem 11.31 

Next, we state Theorems l4.2l and l4.3l from which Theorems ll.5l[L6l and Corollary 
O follow. These results concern a wider range of parameters, and enable us to 
understand something about the border of the basin of attraction of the null solution 
and of infinity, i.e. blowing up solutions. We start from a generalization of a result 
by Wang in m concerning slow decay solutions. 

4.2. Theorem. Assume Fu, Fs, Gu, Gs. Assume either A~^ withlu,ls S (2*, 2*] 
or A~ with lu > 2* and 2* <4 < a*. Then there exists a one-parameter family 
of upper radial solutions of denoted by Xt{x), such that u(x,t]x) converge to 
the null solution as t ^ +00 and with the properties described in Theorem \1.5[ in 
particular they have slow decay. 

Such a result is proved in m for f{u,r) = but as far as we are aware 

is new even for the potential considered in [I] and [55]; be. even for / as in 
with k{r) ^ r^ and / as in (11.41) also when fci = /c 2 = 1- Theorem 14.21 seems to 
suggest that slow decay is the optimal decay rate to have solutions of CH) defined 
for any t. In fact we have the following result which is in contradiction with this 
idea (as far as we are aware these results are new even for the case f{u) = 

4.3. Theorem. Assume Fu, Fs, Gu, Gs. Assume either A'^ with lu, Is G (2*, 2*] 

or A~ with lu,ls > 2*. Then there are one-parameter families of upper and lower 
radial solution of denoted by C,t {x) and ipr (x) having the properties described 

in Theorem [HB 

We stress that Remark 11.71 holds also in this case. 

The relevance of Theorems 14.21 and 14.31 is more clear if we recall the comparison 
principle: if we choose any (j){x) even non radial, such that p (Cr for some r then 
u{x, t] p) converge to the null solution, while if p 'Pt{x) then u(x, t; p) blows up 
in finite time. In fact Corollary II.81 holds in this more general context too. 

4.1. Construction of upper and lo-wer-solutions for the stationary prob¬ 
lem. From now to the end we always assume FO, Fu, Fs (in order to guarantee 
local existence of solutions) without further mentioning. We introduce the follow¬ 
ing notation; we denote by y“(s, a; /„) the trajectory of (12.211 corresponding to the 
regular solution U{r,a) of (11.811 . by y“(s,oo;lu) the trajectory corresponding to 
the singular solution U{r,oo), by y®(s,/3;/s) the trajectory corresponding to the 
fast decay solution V (r, /?), and by y^{s, 00 ; Is) the trajectory corresponding to the 
slow decay solution t7(r, 00 ). 

4.4. Lemma. Fix S G K. Assume lu > 2*, theny'^{s,a;lu) converges to y^(s, 00 ; lu) 
as a ^ -boo, uniformly for s < S. 

Assume 2* < < 2*, then y^{s,l3-,ls) converges to y®(s, oo;/s), uniformly for 

s > S. 
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We think it is worth recalling that in the cases considered U{r,oo) and 1^(7-, oo) 
are the unique singular and slow decay solutions of m- 

Proof of Lemma Assume lu > 2*, fix r € R and set Q{a) = y“(r, We 

recall that W“(r; lu) is a path joining the origin and R = +oo; lu); moreover 
Q{a) —>■ i? as a —?► oo, see Remark [2.61 and [5]. If / C R is a compact interval, 
then y'^{t, a; lu) —t y“('r, +oo; lu) for any t € I. But a; lu) and y'^{T, +oo; lu) 
are solutions of (12.211 hence y'^{t,a;lu) is equibounded and equicontinuous for a 
large so we conclude using Ascoli theorem. The case 2, < 4 < 2* is completely 
analogous. □ 

Analogously, using the fact that y“(s, a; lu) (0, 0) as s ——oo, and y^{s, /?; Ig) 
(0, 0) as s —>■ +0O for any a > 0 /3 > 0, we get the following 

4.5. Lemma. Assume Gu with lu > 2», and fix S € R; then the trajectory 
y^(s, 0 . 2 ; lu) converges to y'^{s, ai; lu) as 0:2 —t ai, uniformly for s < S. 

Assume Gs with Ig > 2», and fix S € R; then r/®(s,/? 2 ;^s) converges to y^{s, I3i;lg) 
as 132 ^ Pi, uniformly for s > S. 

From Lemmas 01 and EM we easily get the following. 

4.6. Lemma. Let p > 0 be arbitrarily small. Assume Gu, Gs with lu, Ig > 2*, and 
. Then for any ai > 0, U(r, a 2 ) converges to U{r, ai) as 02 —t ai, uniformly for 

r > 0. Further, assume lu > 2*, then U[r,a 2 ) converges to U{r,oo) as a 2 —t + 00 , 
uniformly for r > p. 

Similarly assume Gu, Gs with lu,lg S (2*, 2*], and A~^. Then for any Pi > 0, 
V(r,P 2 ) converges to V(r,Pi) as P 2 —t Pi, uniformly for r > p. Moreover, if 
2t. < Ig < 2* then we can take also Pi = + 00 . 

Proof. Assume Gu, Gs with lu,lg > 2*, and A'~ and choose ai G (0,+ 00 ). Then 
from Lemma EM we easily see that C/(r, 0 : 2 ) converges to U(r, ai) uniformly for 
r in compact subsets of r G (0,c»). However using Ascoli theorem and working 
directly on (11.81) we easily get uniform convergence for r G [0, p) too, see e.g. the 
Appendix of [5] for more details (in a much more general context); so we have 
uniform convergence in [0,i?]. From Proposition 12.12] we know that U{r,a) is a 
G.S. with slow decay for any a > 0, hence yf{s,a;ls) is positive and bounded for 
s > 0. Thus setting r = e'* > 1 we find 

|17(r,a2) - U{r,ai)\ = \y^{s,a 2 ;lg) - y^{s,ai;lg)\e-^^^‘^^ < iFr—. 

So for any e > 0 we can choose Rq = [K/8)^/'^^^“'^ so that \Lf{r, a 2 ) — U{r,ai)\ < e 
for r > Rq. Then we can choose R > Rq and we have that t/(r, 02 ) converges to 
U{r,ai) uniformly for r > 0. 

When lu > 2* we simply repeat the argument for U{r,oo) using Lemma 14.41 
instead of Lemma EM 

Now, assume lu,ls G (2,, 2*] and that A+ holds, so that V{r,Pi) is a S.G.S. with 
fast decay. Then we get uniform convergence for r > p directly from Lemma 14.51 
and if Ig <2*, we conclude by using Lemma HU] □ 

4.7. Lemma. Assume A~,Gu,Gs with lu > 2* and Ig G [2*, a*), then for any 
0 < ai < 0:2 < 00 there is Z(a 2 ,ai) > 0 such that Lfir, a 2 ) > U{r, ai) for 
0 <r < Z{a 2 ,ai) and U{r,a 2 ) = U{r,ai), U'{r,a 2 ) < U'{r,ai) forr = Z{a 2 ,ai). 
Assume A+,Gu,Gs with at, < lu < 2* and Ig G (2*, 2*], then for any 0 < Pi < 
P 2 < 00 there is W{P 2 ,Pi) > 0 such that V{r,P 2 ) > V{r,Pi) for r > W{P 2 ,Pi) 
andV{r,P 2 ) = V{r,Pi), V'(r,P 2 ) > V'{r,Pi) for r = W{P 2 ,Pi). 
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Proof. Assume A~, Gu, Gs with lu > 2* and Is G [2*, ct*). Then from Proposition 
12.121 we know that U{r,a) is a G.S. with s.d. and that it is a S.G.S. with s.d. for 
a = oo. Observe that U{r, 02 ) > U{r, ai) for r in a right neighborhood of 0, since 
[ 7 ( 0 , 02 ) = 02 > Q;i = [7(0, oi), and they are continuous functions in r. Since they 
are slow decay solutions, from Remark |2.10l we see that there is i? > 0 (depending 
on oi, 02 ) such that [7(i?,02) = [7(i?,oi). Then we denote by 

(4.1) Z(o 2 ,oi) := min{i? > 0 | U{R,a 2 ) = U{R,ai)} 

Thus by construction we get ■^U{R,a 2 ) < ^[7(i?,oi) for R = Z(a 2 ,oi), but 
from the uniqueness of the solution of Cauchy problem for ODEs we see that the 
inequality is actually strict. 

The case A~^,Gu,Gs with < lu < 2* and 4 G (2»,2*] is completely analo¬ 
gous, and its proof can be obtained repeating the argument for fast decay solutions 
V{r, (3) and reversing the direction of s. □ 


4.2. Proof of Theorem 14.11 In this subsection we assume the hypotheses of 
Theorem 14.II without further mentioning. Let us set 


(4.2) 


f [7(|x|,a2) \i\x\< Z{a 2 ,ai) 

\ u{\x\,ai) if jxj > Z(Q; 2 ,ai) 

f U{\x\,a2) if |a;| > ^(a2,ai) 
\ U(\x\,ai) a \x\ < Z(a 2 ,ai) 


where Z{a 2 ,ai) is dehned in (14.11) . Then by construction C(a;) and 'ip{x) are radial 
Cs-mild upper and lower solutions for (11.71) . From Theorem Id. 101 and Remark id.lll 
it follows that u{x, t; ^) is radial, decreasing in t and converges uniformly to a radial 
non-negative solution u{x, T^; C) of (11.71) for t < Tq. Thus ||M(a;, t; C)(l -|- |a;|'"(*“^)||oo 
is bounded for t < T(^, hence = -|-oo. Since the null solution is the unique 
radial solution of Ha staying below C{x) for any x G K", see Lemma [2.151 then 
limt_>oo ||u(a;,t;C)(l + = 0 for any 0 < ly < m{ls). 

Now let (j) G Cb such that there is 02 > 0 and (j){x) ^ U{\x\,a 2 ). From strong 
maximum principle for parabolic equations (see, e.g. the appendix in [H]), we get 
u{x,t-,<f>) < U{\x\,a 2 ) for any x and any t > 0. So, up to a time translation, we 
can assume (j){x) < U{\x\,a 2 ) for any x G M". From Lemma [4l6l we see that for any 
e > 0 we can hnd oi < 02 such that |[7(|a;|,ai) — U{\x\,a 2 )\ < £ for any |a;| G K". 
Let C(;c) be the upper solution defined by (14.2p . then if £ > 0 is small enough we can 
assume (j){x) < ({x) for any a; G M”. Hence 0 < u{x, t; (j)) < u{x, 7; () for any a; G M" 
and any t > 0; so = -|-oo and ||n(a;, t; <)')||oo —>• 0, ||u(a;, 7; (j)){l + |a;|^)||oo —>■ 0 as 
7 —>■ 00 , for any 0 < < m{ls). 

Similarly consider (f) G Cb and assume that there is oi > 0 such that (j){x) > 
[7(|a;|, oi), and ^(a;) ^ [7(|x|,ai). Reasoning as above we can assume (/)(x) > 
U{\x\,ai) for any x G M", and we can find 02 > cti, with 02 — oti small enough 
so that the lower solution if(x) dehned in (14.21) satishes '0(x) < 4>{x) for any x G 
K”. From Theorem Id.lHI and Remark Id.Ill it follows that u{x,t;'tp) is radial and 
increasing in 7 and converges uniformly to a radial non-negative solution U{x) of 
(II3 for t <T.ip and U{x) > '(/'(a;) for any x. But from Lemma r2.15l we see that such 
a solution [7(a;) does not exist, hence T .0 < 00 and limt-^T,/, \\u{x,t]tp)\\oo =-l-oo. □ 


4.d. Proof of Theorems 14.21 and 14.31 The construction of the family of upper 
and lower solutions of Theorems 14.21 and 14. dl is based on Remark 12.71 Lemmas l2.131 

[HI 

Proof of Theorem Assume hrst A~,Gu,Gs, with lu > 2* and C > 2*. We 
recall that the critical point {P~°°,Q) of (12.101) admits a 1-dimensional unstable 
manifold and that we denote by (y'“(s; 7^), 2 ( 5 )) the unique trajectory belonging 
to this manifold, and by [7(r, 00 ) the corresponding solution of (11.81) which is a 
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S.G.S. with slow decay, and by y“(s; Is) the corresponding trajectory of (12.21) with 
I = Is- From Lemma [2.131 we know that for any r e R, LF“(r;/s) and y’^ir^ls) 
are contained in i.e. the bounded set enclosed by LF®(r;/s) and the 2/2 

coordinate axis (see the construction just before Lemma [2.1,311 . 

Observe that y hence the values 

R{lu) := min{ 2 /“(s; lu) \ s < 0} R{ls) := min{ 2 /“(s; Is) \ s > 0} 

are strictly positive and bounded. Assume first that < Ig, so that m(/„) > m{ls), 
and set 


(4.3) R := 1/2 min ({i?(lti)exp[(m(ls) - m(/„))T] | r < 0} U {R{ls)}) 

Note that R > 0 and y/(T; Ig) > R for any s £ M. 

If lu > 2* then LF“(r; Ig) is a path connecting the origin and y“(T; Ig) so we can find 
(at least) a point Q'^’*{t) = {Qi’*{t),Q 2 ’*{t)) G 1F“(t;1s) such that Qi’*{t) = R. 
Moreover there are two points, say (5®’"*"(t), (5®’“(t), belonging to LF®(t;1s) and 
such that Q®’^(t) = (.R, (^^’^(r)), and Q^ij) < P^i''') < Qt (''')■ Let us consider 
now the trajectories y(s, t; (5®’+(t); /g) and 2 /(s, r; Q®’“(r);/g): by construction 
they correspond to fast decay solutions of (11.81) . say V{r,/32) and V{r,j3i) respec¬ 
tively. We can assume w.l.o.g. that /32 > /?i, see Lemma 1^31 and Remark l^jHl We 
denote by U{r, a*) the regular solution of (11.811 corresponding to y(s, r; Q‘^’*{t)-Is). 
Since 1 / 2 ( 5 , r; (r); Ig) < y 2 (s, r; Q“’*(t); Ig) < y 2 (s, r; Q®’+(t); /g) we have 
V'{R,l3i) < U'{R,a*) < V'{R,fi 2 ) for R = . Therefore we can construct up¬ 

per and lower radial solution of ini) , say C,{x) and ^^{x) as follows: 


(4.4) 


f U{\x\,a*) if |a:| < e"^ 

I y\\x\,lii) if ia:i > e^ 

f 17(|x|, a*) if \x\ > 

\ F(|x|,/ 32) if Ixj < e^ 


Moreover observe that Pi —>■ -boo and a* —)► 0 as r —-boo, and Pi ^ 0 and a* 
-boo as r —— 00 , for i = 1, 2, see RemarkHence C(0, r) = '0(0, r) := D{t) —/■ 0 
as T —>■ -boo and Dir) —>■ -boo as r —>■ — 00 . Similarly lim| 3 ,|_^+oo (/(x, t)|x|"“^ = Pi 
and lim| 3 ,|_>,+oo ■0(x, t) = /32 go to 0 as r —>• -boo and they go to -boo as r —/• 0. 

So, from Theorem 13.101 Remark 13.111 Lemma 12.151 we see that for any r £ R. 
||u(x, t; C(x))||oo —>■ 0 and ||m(x, t; C(x))(l -b |x0)||oo —0 as t —>■ -boo for any v £ 
[0,n — 2), and that there is such that ||u(x, t;'0(x))||oo — +00 as t ^ T^. 

Now we go back to the case where 2* < Ig < lu so that m{lg) > m{lu), and A~ 
holds. In this case we need to consider 


(4.5) R := 1/2 min ({i?(lg)exp[(m(lu) - m{lg))T] | r > 0} U {R(/„)}) 

and to redefine a set E^{t) := {Q = | ^ ^ (note that we 

could redefine E^{t) simply considering the bounded set enclosed by LF®(t; /„) and 
the 2/2 axis), and observe that LF“(1 u;t) C E^{t) for any r £ R. Then we repeat 
the previous argument working with (EH) with I = lu and replacing R by R. 

Now assume A~, with lu = 2* and Ig > 2*. In this case W^ir-Ps) is a 
manifold departing from the origin and contained in E^(t) for any r £ R, but it 
does not have y“(r;Zg) in its border. However for any r £ R we can still find 
at least an intersection between the line yi = R and LF“(T;/g) and at least two 
intersection between yi = R and W‘‘{t-,Is). So it is easy to check that the whole 
argument can be repeated with no changes. 

Now assume A~, with lu = Ig = 2*. In this case we rely on Lemma [2.141 let 
U{r,a*) be the regular solution of (11.81) corresponding to y(s, r, (5“’+(r); 2*), and 
V{r,P 2 ), U{r,Pi) be the fast decay solutions of (11.81) corresponding respectively to 
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y(s, T, (5®’+(r); 2*) and y{s,T,Q^’ (t);2*). Then we define and ipix) as in 
(EH): they are respectively radial upper and lower solutions for (113. Moreover 
C(0,t) = tp{0,T) := D{t) —0 as t —+oo and D{t) —+oo as r —>• —oo. Similarly 
lim| 2 :|_>+oo C{x,t)\x\^~'^ = /3i and lim|a;|_>,+oo 4’ix,T)\x\^~^ = /32 go to 0 as r +oo 
and they go to +oo as r —>■ 0. 

The case where A~^, holds with 2* < < 2* and 2, < Zg < 2* or lu = Is = 2* 

are completely analogous and are left to the reader, so the proofs of Theorem 14.dl 
and Proposition 1 1.71 is concluded. □ 

Proof of Theorem 14-31 Assume A~, then Theorem 14.21 follows from Theorem 14.11 
So assume A+ and assume first Z„,Zs € (2*, 2*). Then from Proposition 12.121 
there is a unique S.G.S. with s.d. say Z7(r,oo) of (11.81) : let y^(s;Z„) be the corre¬ 
sponding trajectory of (12.21) . Fix t £ K: from Lemma 12.131 we know that there is 
Q(t) = (Qi(t),Q 2 (t)) G w,“(r) such that Qi(t) = yf(T;Q and Q 2 (t} < y|(T;Z„). 
Consider the trajectory y(s,T,Q(T);lu) and the corresponding regular solution 
U(r,d(T)) of p.8l) : from Remark 12.61 we see that cZ(r) —>• 0 as r —^ -l-oo and 
cZ(r) — >■ -1-00 as T — >■ —oo. Moreover the function x ( x , t ) defined as follows 


(4.6) 


/ U{\x\,d{T)) if |a;| < e^" 
( Z7(|a;|, oo) if |a;j > 


is a super-solution of dni), and it is regular and has slow decay. Hence, rea¬ 
soning as above we see that u{x, t, x) is radial, radially decreasing and well de¬ 
fined for any t and converges monotonically to the null solution as Z —> -l-oo and 
limt^+oo ||M(a:,t,x)[l -|- |a:|'"]||oo = 0 for any 0 < i/ < m(Zs). If 2 * < Z„ < Zg = 2 *, 
we just lose uniqueness of the S.G.S. with s.d. but the argument can be repeated 
for any such solution, so it still works. If G = G = 2* we simply need to repeat the 
same argument but using Lemma l2 . 141 instead of Lemma [ 2.131 □ 
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